Deconfinement and color superconductivity in cold neutron stars 
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We study the deconfinement transition of hadronic matter into quark matter in neutron star 
conditions in the light of color superconductivity. Deconfinement is considered to be a first order 
phase transition that conserves color and flavor. It gives a short-lived (r ~ T^eak) transitory 
colorless-quark-phase that is not in /3-equilibrium. We deduce the equations governing deconfinement 
when quark pairing is allowed and find the regions of the parameter space (pairing gap A versus 
bag constant B) where deconfinement is possible inside cold neutron stars. We show that for a wide 
region of {B, A) a pairing pattern is reachable within a strong interaction timescale, and the resulting 
"2SC-like" phase is preferred energetically to the unpaired phase. We also show that although /3- 
stable hybrid star configurations are known to be possible for a wide region of the (_B, A)-space, 
many of these configurations could not form in practice because deconfinement is forbidden, i.e. the 
here studied non-/3-stable intermediate state cannot be reached. 



I. INTRODUCTION 

A general feature of degenerate Fermi systems is that 
they become unstable if there exist any attractive inter- 
action at the Fermi surface. As recognized by Bardeen, 
Cooper and Schrieffer (BCS) (Ij this instability leads to 
the formation of a condensate of Cooper pairs and the 
appearance of superconductivity. In QCD any attractive 
quark-quark interaction will lead to pairing and color su- 
perconductivity, a subject already addressed in the late 
1970s and early 1980s 0,0 which came back a few years 
ago since the realization that the typical superconducting 
gaps in quark matter may be larger than those predicted 
in these early works (A as high as ~ 100 MeV) P|. The 
phase diagram of QCD has been analyzed in the light of 
color superconductivity and model calculations suggest 
that the phase structure is very rich at high densities. 
Depending on the number of flavors, the quark masses, 
the interaction channels, and other variables many possi- 
ble /3-stable color superconducting phases of quark mat- 
ter are possible UlElllSli- 

There is at present some indication that the quark 
gluon plasma might have been produced in laboratory 
[lO|. However, it is not yet established whether decon- 
finement happens in nature in the high density, low tem- 
perature regime that is relevant for neutron stars. Un- 
fortunately, first principle calculations are not available 
in this region of the QCD phase diagram. In turn we 
shall base our analysis on phenomenological considera- 
tions which could delineate at least a broad brush pic- 
ture of the physics involved. Matter in compact stars 
should be electrically neutral and colorless in bulk. Also, 
any equilibrium configuration of such matter should re- 
main in /3-equilibrium. Satisfying these requirements im- 
pose nontrivial relations between the chemical potentials 
of different quarks. Moreover, such relations substan- 



tially influence the pairing dynamics between quarks, 
for instance, by suppressing some color superconducting 
phases and by favoring others [TH . 

At not very large densities, the appearance of strange 
quarks is suppressed because of their flnite mass and 
quark matter is composed almost completely by quarks 
u and d. Pairing between them would be possible if 
their Fermi momenta are not very different. The re- 
sulting pairing pattern is the so-called two-flavor color 
superconductor (2SC) in which the up and down quarks 
form Cooper pairs in the color-antitriplct, flavor-singlet, 
spin- zero channel. In the conventional picture of the 2SC 
phase it is assumed that pairs are formed by up-red (u^) 
and down-green {dg) quarks, as well as by up-green {ug) 
and down-red {dr) quarks. The other two quarks {ub and 
dh) do not participate in pairing. A more recent analy- 
sis shows that the ground state of dense up-and-down 
quark matter under local and global charge neutrality 
conditions with /^-equilibrium has at least four possibili- 
ties: normal, regular 2SC, gapless 2SC phases, and mixed 
hase composed of 2SC phase and normal components 
A new interesting feature of some of these phases is 
that pairing is allowed between particles even in the case 
> A 0. 

At sufficiently large densities the value of the chemi- 
cal potential exceeds the mass of the strange quark m^. 
Strange quarks appear in the mixture, and Cooper pair- 
ing can happen between up, down and strange quarks. 
Pairing involving strange quarks is expected to exist if the 



resulting gaps (A„s and A^s) are larger than 
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the difference between the u and s Fermi momenta in 
the absence of pairing. Depending on the value of the 
strange quark mass, as well as other parameters in the 
theory, many different paired configurations are possible. 
At very high densities it seems clear that the color-flavor 
locked phase (CFL) is the ground state. However, at 
not very large densities, it is possible that up and down 
quarks form 2SC matter, while the strange quarks do 
not participate in pairing, eventually forming a < ss > 
condensate with a much smaller gap. Finally, even more 
exotic crystalline phases of ISC quark matter have been 



2 



analyzed. 

To the best of our knowledge, all previous works about 
color superconductivity in compact stars have dealt with 
matter in /^-equilibrium. This is the situation expected to 
appear in strange stars or hybrid stars as soon as they set- 
tle in a stable configuration. However, during the decon- 
fincment transition in neutron stars, matter is transito- 
rily out of equilibrium with respect to weak interactions. 
In fact, the transition from /3-stable hadron matter to 
quark matter in cold neutron stars should occur tr ough 
a quantum nucleation process ilS, iHk dfib lO; ^M- 
Quantum fluctuations could form both virtual drops of 
unpaired quark matter (hereafter the Q*„p phase) or vir- 
tual drops of color-superconducting quark matter {Q*^ 
phase). In both cases, the flavor content of the quark 
matter virtual drop must be equal to that of the confined 
/3-stable hadronic phase at the same pressure (the central 
pressure of the hadronic star). In fact, since quark decon- 
finement and quark-quark pairing are due to the strong 
interaction, the oscillation time of a virtual quark 
droplet in the potential energy barrier separating the 
hadronic from the quark phase, is of the same order of the 
strong interaction characteristic time (Tstrong ^ 10~^^ s). 
The latter is many orders of magnitude smaller than the 
weak interaction characteristic time {Tweak ~ 10~^ s). 
Thus, quark flavor must be conserved forming a virtual 
drop of quark matter d El El IM 113 • 

Which one of 

the two kind of droplets {Qunp Q*a) will nucleate de- 
pends on the value of the corresponding Gibbs free energy 
per baryon {gunp, 5a)- In fact, the latter quantity enters 
in the expression of the volume term of the energy barrier 
separating the confined and deconfined phases (see e.g. 
eq. (7) in 18], where the Gibbs free energy per baryon is 
denoted by fii^i — Q*,H). Clearly, when < gunp the 
nucleation of a Q\ drop will be realized. 

The direct formation by quantum fluctuations of a 
drop of /3-stable quark matter (Q phase) is also possi- 
ble in principle. However, it is strongly suppressed with 
respect to the formation of the non /3-stable drop by a 
factor ~ Gp^{^; being TV the number of particles in the 
critical size quark drop. This is so because the formation 
of a /5-stable drop will imply the almost simultaneous 
conversion of ^ iV/3 up and down quarks into strange 
quarks. For a critical size /3-stable nugget at the center 
of a neutron star it is found N ^ 100 — 1000, and there- 
fore the factor is actually tiny. This is the same reason 
that impedes that an iron nucleus converts into a drop of 
strange quark matter, even in the case in which strange 
quark matter had a lower energy per baryon (Bodmer- 
Witten-Terazawa hypothesis). Because of this reason it 
is assumed that a direct transition to /3-stable quark mat- 
ter is not possible [3l|. Therefore, the /3-stable state Q 
could be reached only after the /3-decay of the interme- 
diate state Q* . This is in agreement with many other 
previous works, see e.g. [Ulil El El EH . 

In this context, one question addressed by the present 
work is whether the system settles in a paired or in an un- 
paired state just after the deconfinement. On the other 



hand, notice that although the above mentioned non-/3- 
stablc quark phase is very short-lived, it constitutes an 
unavoidable intermediate state that must be reached be- 
fore arriving to the final /3-stable configuration, e.g. CFL 
quark matter (c.f. (22,,^J). The second question we shall 
address is whether this intermediate phase can eventually 
preclude the transition to the final /3-stable state in spite 
of the latter having a lower energy. This is because the 
Q-phase can be formed only after the nucleation of a real 
(i.e critical size) drop of Q*„p or Q\ matter, and its sub- 
sequent "long term" {t ~ T^eak lO^^s) weak decay 
process. 



II. DECONFINEMENT OF HADRONIC 
MATTER INTO COLOR SUPERCONDUCTING 
QUARK MATTER 

Given the uncertainties in the nature of matter at high 
densities, the analysis is based on the extrapolation to 
higher densities of an hadronic model valid around the 
nuclear saturation density poi f^nd the extrapolation to 
Po of a quark model that is expected to be valid only for 
p — > oo. Within this kind of analysis the (in general) dif- 
ferent functional form of both EOSs, induces the phase 
transition to be first order. Notice that from lattice QCD 
calculations there are indications that the transition is ac- 
tually first order in the high-density and low-temperature 
regime, although this calculations involve temperatures 
that are still larger than those in neutron stars, and do 
not include the effect of color superconductivity . 

Deconfinement is analyzed here as a first order phase 
transition that conserves the fiavor abundances in both 
phases. Therefore, it gives a transitory colorless-quark- 
phase that is not in /3-equilibrium. 

For describing the just deconfined quark phase we shall 
model it as a free Fermi mixture of quarks and leptons 
and we will subtract the pairing and the vacuum energy. 
The thermodynamic potential can be written as 

VI = r!^^^ + + f^L + B, (1) 

with 



"q" = E4 / \E,f - Nc]p^dp, (2) 

where Ecf = {p^ + m^^)^/^, and the sum is to be made 
over all colors and fiavors of the quark mixture, being 
/ = u,d,s the flavor index, and c — r,g,b the color 
index. Confinement is introduced in Eq. |^ by means 
of the bag constant S, and fii is the contribution of 
leptons. Pairing is included through the term 

= (3) 
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FIG. 1: Schematic comparison of the free energy of hadronic 
matter (H), non-/3-stable "just-deconfined" matter (Q*), and 
/3-stable quark matter (Q) for different cases. In the case of 
panel a) the transition can never occur inside neutron stars in 
spite of the final state (Q) having a lower energy per baryon. 
As explained in the text, a direct transition to Q is strongly 
suppressed. Since Q* has a larger energy per baryon than H 
for all pressures below the central pressure of the maximum 
mass hadronic star (P™''"), deconfinement cannot occur even 
if the Q phase has a lower energy. In panels b) and c) the 
phase Q* has a lower energy per baryon than H for some 
pressures below (P™'"^). Therefore, deconfinement is possible 
if pressures between Pq and P™''" are reached inside a given 
neutron star. The difference between panels b) and c) is the 
energy per particle at zero pressure (indicated with a dot for 
Q and with an asterisk for H). In b) quark stars are the so 
called strange stars, because they can be made up of quark 
matter from the center up to the surface (P=0). In case c) 
they are hybrid stars, because at zero pressure H has a lower 
energy than Q. For a fixed hadronic equation of state, these 
possibilities correspond to different values of the parameters 
of the quark model (the vacuum energy density B and the 
superconducting gap A). 



which results from Ecf — \/p'^ + m"^ ± A^, by expand- 
ing in Eq. ((2Jl to 0{A'^/^'^). For simpHcity, all particles 
that pair are assumed to have the same gap A. The 
mean chemical potential in Eq. |(2Jl is defined through 
p, = c /^/c where the sum is to be made only over 

particles that participate of pairing, and N is the num- 
ber of different quarks that pair. The Fermi momenta 
are kfc = {iJ.^f^ — m^^)^/^. The binding energy of the di- 



for every quasiparticle with gap A [23. As we shall see 
below, the relevant pairing pattern for the deconfinement 
transition in neutron stars is "2SC-like" , in which dr and 
Ug quarks pair yielding two quasiparticles with gap A, 
and Ur and dg quarks pair, yielding two quasiparticles 
with gap A (c.f. [l^). Therefore, we shall set ^ = 1 
along this work. Note that for color-flavor-locked quark 
matter it is „4 = 3. 

As already emphasized by Rajagopal and Wilczeck 0, 
the exact nature of the interaction that generates A is 
not relevant to the order we are working. This means 
that ri is given by the above prescription regardless of 
whether the pairing is due to a point like four-Fermi in- 
teraction, as in Nambu-Jona-Lasinio models or due 
to the exchange of a gluon, as in QCD at asymptoti- 
cally high energies. Of course, the strength and form of 
the interaction determine the value of A, and also its 
dependence with the density. Lacking of an accurate cal- 
culation for A, which may as large as ^ 100 MeV (and 
even larger in the presence of an external magnetic field 
), we shall keep it as a free constant parameter. 

The thermodynamic quantities are straightforwardly 
derived from the standard expressions: the pressure is 
P = —fl and the energy density at zero temperature is 
given by £ = ^ /x^ni -I- O, where the sum is to be 
carried over all quarks and leptons. In particular, the 
particle number densities nfc = —dft/d^fc are given by 
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2A 



37r2 Ntt^ 

for the quarks that participate in pairing, and by 



3^ 



(4) 



(5) 



for quarks that do not pair. The number density of each 
flavor in the quark phase is given by 



and the baryon number density ti^ by 



(6) 



(7) 



fc 



quark condensate is included by subtracting A 11 / {'in ) 



In order to close the above equations we need to impose 
additional physical conditions describing the composition 
of the mixture (e.g. a set of conditions on the chemical 
potentials). One possibility, extensively employed in the 
literature, is /?-equilibrium of the quark phase. This con- 
dition describes matter a sufRciently large time after de- 
conflnement (r 3> Tweak ~ 10~^s). However, for study- 
ing deconfinement, the relevant timescale is r <C T^eak- 
As already emphasized in IHIll [H IIH the ap- 
propriate condition for r <C T^eak is flavor conservation 
between hadronic and deconfined quark matter. This can 
be written as 
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yH 



f — u^d, s, L 



(8) 



being Yy^ = n^^ /n|| and Y^*^ = nj /n^ the abundances 
of each particle in the hadron and quark phase respec- 
tively (we shall omit the super-indexes H and Q in the 
following). In other words, the just deconfined quark 
phase must have the same "flavor" composition than the 
/3-stable hadronic phase from which it originated. 

Additionally, the deconfined phase must be locally col- 
orless; therefore, it must be composed by an equal num- 
ber or red, green and blue quarks: 
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llr — Hg — ni, (9) 

being n^, Ug and rib the number densities of red, green 
and blue quarks respectively, given by: 

ric = J^^f'^- (10) 
/ 

Color neutrality can be automatically fulfilled by im- 
posing that each flavor must be colorless separately, i.e. 

^ur — ^ug — ^uh} ^dr — ^dg — ^dbj and TL^y — ^sg — ^sh- 

But in general this configuration will not allow pairing 
with a significant gap. As already stated, for quarks hav- 
ing different color and flavor the pairing gap may be as 
large as 100 MeV, while for particles having the same 
flavor the gap is found to be about two orders of magni- 
tude smaller (see and references therein). Pairing is 
allowed even in the case (5/i > A but the corresponding 
gaps are small [^ll^. Therefore, in order to allow pair- 
ing between quarks with a non negligible gap, the Fermi 
momenta of at least Ur and dg quarks must be equal (the 
choice of these two particular colors and flavors is just a 
convention). This implies the equality of the correspond- 
ing number densities 

nur = ridg- (11) 

The above condition represents a state that fulfills all the 
physical requirements of the deconfined phase (e.g. is 
color and electrically neutral), and should be the actual 
state (for r <C Tweak) if it has the lowest free energy per 
baryon. Energy must be paid in order to equal at least 
two Fermi seas, but in compensation the pairing energy 
is recovered. The gained energy depends on the value of 
the pairing gap A and, at least for sufficiently large A, 
is expected to be larger than the energy invested to force 
a pairing pattern. Also, notice that color conversion of 
quarks allows the adjustment of the Fermi seas within 
a given flavor in a very short timescale (~ Tstrong), i-e. 
several orders of magnitude faster than /3-equilibration 

{^strong Tweak)- 

We emphasize that this phase is not in flavor equilib- 
rium. After a weak interaction timescale this transitory 
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FIG. 2: Deconfinement of pure neutron matter. Upper panel: 
Sketch of the lowest energy configuration of the paired and un- 
paired phases just after the deconfinement. Lower panel: The 
difference in the Gibbs energy per particle between paired and 
unpaired quark matter. For positive values of qa — Qunp the 
preferred phase just after the deconfinement is the unpaired 
one while for negative values it is the paired one. 



pairing pattern will be abandoned by the system in favor 
of the lowest-energy /3-stable configuration. Depending 
on the density, the lowest energy state may be LOFF, 
gapless 2SC, gapless CFL, standard CFL, (to name just 
some possibilities) as extensively discussed in the litera- 
ture. 



III. APPLICATION TO SIMPLIFIED 
EQUATIONS OF STATE 

A. Deconfinement of pure neutron matter 

A simple solution can be found in the case of the de- 
confinement of pure neutron matter, since strange quarks 
and electrons are not present in the hadronic gas. First, 
we apply the colorless conditions and flavor conserva- 
tion introduced in the previous section, in order deter- 
mine the abundances of each quark species. In order 
to allow pairing of at least two different quarks species 
in the just deconflned phase, we impose the condition 
of Eq. (|ll|l . i.e. n^r = n^g- Using one of the color- 
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less conditions (n^ = Ug) it is found that Udr = n^g, 
implying that these quarks can also pair with a sig- 
nificant gap. From the remaining colorless condition 
[rir — rib) it is found + i^dr = nub + n-db- The con- 
dition of flavor conservation states Ud = 2nu', therefore, 
nd7- + ridg + Udb = 2(n„.r + Uug + Uub)- Introducing the 
ratio X — Uug/riur we find from the above equations: 



riub = 



ndb 



1 



(12) 



(13) 



Therefore, for massless particles at T = the chemical 
potentials are related by: 



fJ-ub ^ 



(14) 

(15) 
(16) 



With this configuration, the pressure Pa and the Gibbs 
free energy per baryon = (X]/c "'/cM/c)/"'S take the 
simple form 



Pa 



[2(1 + :e4/3) + (1 + ^)4/3]^. 



4 

ur 



9A 



127r2 
2^^ ^' 



[2(l-Ha;4/3) + (l + x)4/3]/ 



(17) 



(18) 



Minimizing (Ja = 3a(Pa, x) with respect to x it is found 
that the minimum correspond to a: = 1. Therefore, 
quarks Ur — dg and Ug — dr pair in a "2SC-like" pattern 
like the one shown in Fig. [21 

In order to determine whether the system settles in a 
paired or in an unpaired configuration we compare the 
Gibbs free energy per baryon of the above configuration 
with the Gibbs free energy per baryon of an unpaired 
quark gas (both evaluated at the same pressure, and with 
the same fiavor composition). For unpaired quark mat- 
ter the ground state of the colorless mixture (compatible 
with flavor conservation) is shown in Fig. |21 and is de- 



scribed by Udr = TLdg 



ndb 



2n,, 



2n„ 



^ub- The Gibbs free energy per baryon reads 

g„„^ = [471^(1 -H24/3)3(p + B)]l/4. 



(19) 



The results are shown in Fig. [3 where the difference 
9 A — gunp is plotted as a function of pressure for different 
values of the parameter A. For large enough values of 
A the energy of the pairing gap is able to compensate 
the increase of the free energy resulting from equating 
the Fermi seas of the particles that pair. It is also worth 
noting that since we are comparing the Gibbs energy per 
baryon of different phases at equal pressures (alterna- 
tively, we could compare the pressure at fixed Gibbs en- 
ergy per baryon), the results depend on the bag constant 
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FIG. 3: Deconfinement of n-p-e~ gas. Upper panel: Sketch of 
the lowest energy configuration of the just deconfined paired 
and unpaired phases . Lower panel: The difi'erence in the 
Gibbs energy per baryon between both phases at ^ = 1.6. 
For positive values of (/a — gunp the preferred phase just af- 
ter the deconfinement is the unpaired one while for nega- 
tive values it is the paired one. The results are shown for 
B = 60 MeV fm"^ (dashed line), B = 100 MeV fm"^ (solid 
line), and B = 140 MeV fm~^ (doted line). 



B. This dependence is not found in the comparisons for 
/3-stable quark matter made in 12^ because the bag con- 
stant B cancels out when comparing AP at fixed /z. Note 
that AP is not the relevant thermodynamic function for 
the analysis of the deconfinement transition in neutron 
stars performed here. 



B. Deconfinement of n-p-e gas 

Next we consider the deconfinement of a charge neu- 
tral uniform system of neutrons, protons and electrons 
{rim and are the associated number densities). As 
discussed above, we impose n^r — ridg in order to allow 
pairing. Then, color charge neutrality imposes ridr — riug 
and riur + ridr — ri^b + ndb- The resulting configuration 
allows pairing of quarks Ur and dg with a gap A, and of 
Ug and dr with a gap that is assumed to have the same 
value. 

On the other hand, we can put the condition of flavor 
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conservation in the following simple form 



rid = 



(20) 



being ^ = /Y^ (c.f Eq (8)). In the case of the n- 
p-e~ gas, the parameter ^ can be expressed in terms of 
the proton fraction Yp = Up/ns of nuclear matter as 
C = (2 - Yp)/{1 + Yp). It is easy to check that ^ = 2 
corresponds to the deconfinement of pure neutron mat- 
ter, ^ = 1 to symmetric nuclear matter and ^ = 0.5 to 
the unrealistic case of pure proton matter. We emphasize 
that in the case of /9-stable n-p-e" system, ^ is a function 
of density (or pressure) that depends only on the state 
of the hadronic matter that deconfines. 

Therefore, flavor conservation states that n^r + ridg + 
^db = S,{nur + n^g + n^b)- In addition, the condition of 
Eq. © applied to electrons yields: 



3ne — 2nu — Ud, 



(21) 



which confirms that flavor conservation automatically 
guarantees electric charge conservation. Finally, we im- 
pose that 1) ridr = nug in order to allow for paring be- 
tween quarks dr with Ug, and 2) Uur — ridg in order to 
allow for paring between quarks u,. with dg. 

Introducing the ratio x = Uug/riur, and using Eqs. 
H2()|l and H21|l we find the particle number densities of 
each flavor and color in the paired phase as a function 
of one of the particle densities (e.g. n„r), the parameter 
^ that depends only of the state of the hadronic phase, 
and the free parameter x: 



riug 


— IT Tlur 


nub 


= (1+.T 


ndr 


— X Tliif 


ndg 


— ^ur 


ndb 


= (1+.T 



2-e 



2i- 1 



(22) 

(23) 

(24) 
(25) 

(26) 



Note that the free parameter x that can be elimi- 
nated by minimizing the Gibbs energy per baryon = 
fc'^^fcP'fc + fJ'ene)/nB with respect to x at constant 
pressure Pa ■ The minimization gives x = 1, which means 
that the configuration of the paired phase that is energet- 
ically preferred is the one having ndr = n^g = nur = ndg, 
as sketched in the upper panel of Fig. [21 In Fig. O we 
compare the Gibbs energy per baryon of the "2SC-like" 
paired phase and unpaired matter for different values of 
the bag constant B and the pairing gap A. Gompar- 
ing with the pure neutron matter case of Fig. |21 it can 
be noticed that an increase in the proton fraction of the 
hadronic phase favors the formation of a paired quark 
phase after deconfinement. 
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FIG. 4: Deconfinement of /3-stable hadronic matter including 
strange hadrons. Upper panel: Sketcli of tlie particle num- 
ber configuration just after deconfinement. The most general 
paired configuration compatible with the condition Uur = Udg , 
riug ~ ndr, flavor conservation, and charge neutrality is the 
one sketched here. Lower panel: The difference in the Gibbs 
energy per baryon between both phases. For positive values of 
Qa — Qunp the preferred phase just after the deconfinement is 
the unpaired one while for negative values it is the paired one. 
The curves correspond to B = 60 MeV fm~^ (dashed line), 
B = 100 MeV fm"^ (solid line), and B = 140 MeV fm"^ 



(doted line). We employed ms 



150 MeV. 



IV. DECONFINEMENT OF COLD HADRONIC 
MATTER 



In the following we analyze the deconfinement of a gen- 
eral hadronic system including strange hadrons and then 
we apply the results to a realistic EOS in order to study 
deconfinement inside cold neutron stars. 



Deconfinement of a general hadronic equation 
of state 



• Flavor conservation: After deconfinement the parti- 
cle densities of quarks u, d and s are the same as in the 
hadronic phase and can be determined by Eqs. ©. An- 
other equivalent way of expressing the fiavor conservation 
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condition is in terms of two parameters ^ and ry: 

nd = £.nu. (27) 



77 n„. 



(28) 



where ^ = Yj^ /Y^ and 77 = Y^ /Y^ depend only on the 
composition of the hadronic phase. These expressions are 
vaHd for any hadronic EOS. For hadronic matter contain- 
ing n, p, A, S+, YP , E~, S~, and S", we have 



2np + n„ + riA + 2nE+ + n^o + n^o ' 

"A + »lS+ + JT-sri + TIE- + 2nH0 + 2ns- 



(29) 
(30) 



As typical values, we notice that 77 = corresponds to 
zero strangeness, and that at the center of the maxi- 
mum mass star (calculated with the hadronic equation 
of state of Glendenning and Moszkowski GMl |31|) we 
have ^ = 1.15 and 77 = 0.85. Notice that and 77 deter- 
mine univocally the number of electrons present in the 
system through electric charge neutrality of the decon- 
fined phase: 



nd 



(31) 



• Pairing condition: As made in the previous section 
for the n-p-e~ gas, we impose that 1) n^r = nug in order 
to allow for paring between quarks dr with Ug, and 2) 
Uur = ndg in order to allow for paring between quarks u^. 
with do- 



Ug leads 



• Color neutrality: The condition 
immediately to Also, Ur — Uh leads to 

2nur + risr = nub + ndt + ngt- 

The above conditions lead to the pairing pattern 
schematically shown in the upper panel of Fig. ^ Note 
that these conditions still leave a degree of freedom that 
can be fixed by introducing an additional parameter h 
relating the particle number densities of two arbitrary 
quark species. Therefore, it is possible to impose the 
equality of two arbitrary Fermi seas in order to allow 
pairing between them. We have analyzed the 10 possible 
combinations and verified that 6 of them lead to a nega- 
tive value of the particle number density of at least one 
quark species. The other 4 possibilities allow pairing of 
particles that don't have different color and flavor, allow- 
ing pairing with a negligible gap. For this reason, it is 
more convenient to introduce h = Usb/nsr, and minimize 
the free energy with respect to h. 

Using the above Eqs. we find the following linear set 
of equations: 

2nur + nsr = nub + ndb + nsb (32) 
2nur + ndb = ii'^nur + nub) (33) 
2nsr + h nsr = ?y(2n„r + f^t^b), (34) 



from which we obtain the number densities of each quark 
species in the paired phase as functions of only four quan- 
tities: 



" 2-rj + h + 7jh + 2^ + h^ ^-^^^ 
^ -2 + r]-h-vh + 4^ + 2h^ 

2 — 77 -I- ft -I- 77/1 + 2^ -I- /if 

nsb = 7: — ;— — , , , , ^ n^r (37) 

(38) 
(39) 
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-r;4 


- h + t]h 


+ 2^ + 


K 






6?7 






2 


-7^4 


- h + r]h 


+ 2C + 






2(2- 


^h){2- 






2 


-774 


- h + i]h 


+ 2e + 





Remind that the other particle densities are given by 
n^g — ndr — ndg n.ur ^ nsg — n^rj and n^b — hn^r- 

The pressure and Gibbs energy per baryon of the 
paired deconfined phase can also be written in terms of 
the same parameters: 



Pa 



9A 



E 

fc 

E 

fc 



127r2 ^ 

ns 



Me 



127r2 



where kfc = {fJj 
tentials /z/c are obtained from 



B, (40) 
(41) 

the chemical po- 



2A^fL 



fc = Ur, Ug, dr, dg{A2) 



M/c = (37r2,^^c)'/' fc = Ub,db (43) 

[(37r27i^,)2/3+™2jl/2 f^^s^^Sg,Sb (44) 

The minimization of with respect to h gives h — 
1 and therefore the number densities are given by the 
following equations: 



nub = 



ndb 



4- 2C 



nsb 



1 + 




4C- 


- 2 


1 + 
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1 + e 

2(2-77-^) 
1+^ 



(45) 
(46) 
(47) 
(48) 



with Hug = ndr = ndg = nur and Ugg = nsr = "-sb- 
Equations H40|l - H48() constitute the equations of state for 
just deconfined quark matter. In the lower panel of Fig. 
^ we show Ag for particular values of the parameters 
(e = l-6, ry = 0.3). 
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FIG. 5: The cparameter space A vs. B indicating the re- 
gions for which deconfinement is possible inside the maximum 
mass neutron star with the GMl EOS 3 {M^ax = 1.8Mq). 
We also indicate whether the final state reached after l3- 
equilibration of the just-deconfined phase has energy per 
baryon less or greater than the neutron mass (i.e. leads to 
the formation of strange stars or hybrid stars respectively). 
We adopted rris = 150 MeV for the strange quark mass. If 
{B, A) fall inside the dashed region, deconfinement is not pos- 
sible even at the center of the maximum-mass star with this 
EOS. For (B, A) inside the grey region the just-deconfined 
unpaired phase has always less energy per baryon than the 
just-deconfined paired phase. For {B, A) in the white region 
the just-deconfined phase is always paired quark matter. The 
regions met at a point of coordinates {B^,, A«) indicated with 
an asterisk and shown in Table I for different values of the 
strange quark mass m^. The maximum of the grey region is 
indicated with a dot, and the corresponding value Amax is 
shown in Table I for different nis- 



rus (MeV) 


B, (MeV fm"^) 


A, (MeV) A,^ax (MeV) 





323 


34 


78 


100 


300 


37 


80 


150 


275 


49 


85 


200 


241 


68 


97 



TABLE L We show the coordinates of the point (_B«,A») 
(indicated with an asterisk in Fig. |3 for different values of 
the strange quark mass rUs . We also list the coordinate A,nax 
of the point indicated with a dot in Fig. |5] The qualitative 
shape of Fig. |3 remains the same in all cases, but the size of 
each region changes according with the representative points 
given here. 



B. Deconfinement transition in neutron stars 

The above conditions allow us to study the regions of 
the parameter space A versus B where the deconfinement 
transition is possible inside neutron stars (see Figs. [Sjand 
O . In this analysis we shall employ the equations of state 
of Glendenning and Moszkowski GMl for the hadronic 



FIG. 6: The same as Fig. I^but for a 1.6 Mq neutron star. 
The solid curve is the same as in Fig. |5]but the dashed curve 
is strongly shifted to the left, i.e., the region of the parameter 
space that allows deconfinement for a 1.6 Mq neutron star is 
much smaller than for a 1.8 Mq neutron star. Notice that the 
larger part of the available parameter space would lead to the 
formation of strange stars rather than hybrid stars. 

phase Depending on the value of B and A there 

are three possibilities: 1) deconfinement is not possible 
at the center of the neutron star, 2) deconfinement to an 
unpaired phase is preferred, and 3) deconfinement to a 
paired phase is preferred. These regions are limited by 
the following curves: 

• A curve along which the three phases have the 
same state, that is, g_f/(P, 77,^) = gunp{P,riiCj = 
5a(-P, ?7,C) evaluated at the same pressure Ph = 
Punp — Pa and with the same "quark composi- 
tion", i.e. r]H = Vunp = VA and = Cunp = S,A- 
This gives the solid line in Figs. [Sland|Slseparating 
the grey and white regions. 

• A curve along which g, P, ^ and rj of quark mat- 
ter has the same value than at the center of the 
hadronic neutron star (in Figs. |S1 and dashed 
curve for paired quark matter and solid vertical line 
for unpaired quark matter). The position of this 
curve depends on the mass of the neutron star. In 
Fig. |31we show the results for the maximum-mass 
hadronic star within the GMl EOS (1.8 Mq), and 
in Fig. 0for a neutron star with 1.6 Mq. 

The regions meet at a point of coordinates (_B»,A,) 
indicated with an asterisk in Figs. |Sland|^ The position 
of this point (that characterizes rather well the size of the 
regions) depends on the assumed value of the strange 
quark mass rUs, and on the mass of the neutron star. 
In Table I we show the dependence of {B^,,A^) on the 
assumed value of the strange quark mass mg for the 1.8 
AIq neutron star of Fig[Sl We also indicate with a dot 
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the maximum of the grey region (which is the same in 
Figs. [Sjandinjl and give the corresponding value l^max hi 
Table 1. 

We have also included in the parameter space the curve 
separating the regions in which /3-stable quark matter has 
an energy per baryon smaller than the neutron mass from 
the region in which e/nsiP — 0) > m„ (for simplicity, 
paired /3-stable quark matter is assumed in all cases to 
be CFL). To the left of this curve the final state after 
/?-equilibration is absolutely stable quark matter leading 
to the formation of strange stars. To the right, /^-stable 
quark matter is restricted to the core of neutron stars 
(hybrid stars). The position of this curve also depends 
on the value of m^. In Figs. [S] and |B1 it is shown for 
nis = 150 MeV (for more details the reader is refereed to 
0). 



V. DISCUSSION 

In this paper we have analyzed the deconfinement tran- 
sition from hadronic matter to quark matter, and investi- 
gated the role of color superconductivity in this process. 
We have deduced the equations governing deconfinement 
when quark pairing is allowed and, employing a realis- 
tic equation of state for hadronic matter, we have found 
the regions of the parameter space B versus A where the 
deconfinement transition is possible inside neutron stars. 
The main results are shown in Figs. |S1 and and were 
explained in the last section. In the following we discuss 
some implications for neutron star structure. 

Stars containing quark phases fall into two main 
classes: hybrid stars (where quark matter is restricted 
to the core) and strange stars (made up completely by 
quark matter) . This structural characteristic depends on 
whether the energy per baryon of /3-equilibrated quark 
matter at zero pressure and zero temperature is less than 
the neutron mass (the so called "absolute stability" con- 
dition). In the absence of pairing, quark matter in /3- 
equilibrium has an energy per baryon (at P — 0) smaller 
than the neutron mass only if B is in the range 57 MeV 
fm-3 < B < 90 MeV fm^3_ -yVithin this range of B, 
unpaired /3-stable quark matter is the so called strange 
quark matter, and it is possible the existence of stars 
made up entirely by the quark phase. For B > 90 MeV 
fm~'^ unpaired /3-stable quark matter at P = and T = 
decays into hadrons, and therefore it can be present only 
in the core of neutron stars. The size of the core (if any) 
depends on the value of B: the larger the value of B, 
the smaller the size of the quark matter core (for a given 
neutron star mass). 

Pairing enlarges substantially the region of the param- 
eter space where /3-stable quark matter has an energy 
per baryon smaller than the neutron mass ■ Al- 

though the gap effect does not dominate the energetics, 
being of the order (A/yit)^ ^ a few percent, the effect is 
substantially large near the zero-pressure point (which 
determines the stability and also the properties of the 



outer layers and surface of the star). As a consequence, a 
"CFL strange matter" is allowed for the same parameters 
that would otherwise produce unbound strange matter 
without pairing '29*1. The line separating strange mat- 
ter from non-absolutely stable quark matter is shown in 
dotted line in Figs. |31and|Sl according to [2^ . 

Concerning just deconfined quark matter (i.e. not in 
/3-equilibrium) it has been already shown that the transi- 
tion to unpaired quark matter is not possible in a 1.6 Mq 
neutron star if the Bag constant is _B > 126 MeV fm~^, 
because the transition pressure is never reached inside 
the star, even in the proto-neutron star phase The 
results when pairing is allowed have been shown in the 
previous section, where we have shown the "deconfine- 
ment" parameter space for the maximum mass neutron 
star with the GMl EOS (1.8 Mq), and for a 1.6 Mq 
neutron star. As it is evident from Figs. |5l and (HJ de- 
confinement is facilitated for large A (i.e. it is possible 
for a larger range of B). This result can be roughly un- 
derstood if we think paired matter as unpaired matter 
with an effective bag constant depending on the chemi- 
cal potential (or on density): i3off(A, ^) = B — ^A^/i^. 
The minus sign of the condensation term in the previ- 
ous expression allows deconfinement for larger values of 
B. Nevertheless, notice that this simple interpretation is 
not strictly correct because the chemical equilibrium is 
different in paired and unpaired phases. In fact, due to 
the more convenient chemical equilibrium, the unpaired 
phase is energetically preferred for small A. 

For a large part of the deconfined parameter-space A 
vs. i?, /3-stable quark matter has (at any pressure) an 
energy per baryon smaller than the neutron mass. The 
situation corresponds to the one sketched in Fig. lb. 
Therefore, if quark stars are formed, they would be made 
up of quark matter from the center up to the surface. In 
Figs. |5l and El this corresponds to the part of the grey 
and the white regions to the left of the dotted line. 

To the right of the dotted line of Figs. |S1 and [HI stars 
containing quark phases are hybrid. The situation corre- 
sponds to the one sketched in Fig. Ic, i.e. the hadronic 
phase is preferred at low pressures. From the point of 
view of the structural properties, stable hybrid stars have 
been found to be possible in a wide region of the parame- 
ter space [S^- Nevertheless, notice that these configura- 
tions could not form in practice if they fall in the region 
where deconfinement is forbidden (dashed region of Figs, 
m and That is, for a given hadronic star, there exist a 
stable hybrid star with the same baryonic mass that has 
a lower energy (gravitational mass). Nevertheless, the 
hadronic star cannot deconfine because the here stud- 
ied non-/3-stable intermediate state has a larger energy 
per baryon, as shown schematically in Fig. la. This is 
the case, for example, if B = 153 MeV fm'^ {B'^/'^ = 185 
MeV) as can be seen by comparing with the results of Fig. 
5 of the paper by Alford and Reddy : for rus = 150 
MeV, B = 153 MeV fm""^ and any reasonable value of 
A, stable hybrid configurations with maximum masses up 
to 1.6 Mq are found in 32]. However, in these cases the 
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pairs {B, A) fall comfortably inside the dashed region of 
Fig. Elwhere deconfinement is not allowed. For the same 
value of B, heavier stars (~ 1.8 Mq) could deconfine, 
since {B, A) would be inside the grey region of Fig. |S1 
but the resulting configuration would be not structurally 
stable and would form a black hole (c.f. 32]). Although 
the EOSs are different in [s^ and in the present work, 
this should not affect this generic trend. Notice that 
qualitatively similar results have been found in 0, |^ 
for unpaired quark matter. 

A stated in the Introduction, the transition from nu- 
clear matter to quark matter proceeds by bubble nucle- 
ation. However, notice that for large B the results with 
typical surface tension = 10 — 30 MeVfm^^ do not dif- 
fer much from the case in bulk r?,^. This means that 
we don't expect that the dashed region of Figs. [S] and 
l^lwill change significantly when including surface effects. 



Anyway, even if the surface tension were very large, the 
here presented bulk case is still relevant because it gives 
a lower limit for the transition: i.e., if deconfinement is 
not possible in bulk, it will be even more difficult when 
including surface effects. In other words, the dashed line 
of Figs. El and El could move to the left in a more re- 
fined study, but not to the right. A complete study of 
the astrophysical implications is in progress and will be 
published elsewhere. 
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